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Abstract 



We discuss the implications of using an intrinsic Hamiltonian in theories without particle-number conservation, e.g., the Hartree- 
Fock-Bogoliubov approximation, where the Hamiltonian's particle-number dependence leads to discrepancies if one naively re- 
places the particle-number operator by its expectation value. We develop a systematic expansion that fixes this problem and leads 
to an a posteriori justification of the widely-used one- plus two-body form of the intrinsic kinetic energy in nuclear self-consistent 
field methods. The expansion's convergence properties as well as its practical applications are discussed for several sample nuclei. 

fCey words: nuclear structure, intrinsic Hamiltonian, Hartree-Fock-Bogoliubov, particle-number conservation 
PACS: 21.60.Jz 



ON 

o 
o 

> 

o 

£: 

{NJ 1. Introduction 



O 
=5 



> 

m 
m 

00 

o 
o 



X 



Symmetry breaking is a powerful concept in quantum many- 
body theory. In nuclear phyiscs, well-known examples are the 
breaking of the rotational and translational symmetries of the 
many-body Hamiltonian by the use of localized single-particle 
states in the construction of the many-body Hilbert space; the 
latter, in particular, can cause sizable center-of-mass contami- 
nations of the energies and the many-body wave functions un- 
less the symmetry is restored (see e.g. [11120 and Refs. therein). 
In contrast, the breaking of the particle- number symmetry 
by quasi-particle methods like the Hartree-Fock-Bogoliubov 
(HFB) Ull approach is a useful tool because it leads to an ef- 
ficient treatment of the important nuclear pairing correlations, 
although one will ultimately want to restore this symmetry in a 
finite system like the nucleus. 

■ Since nuclei are self-bound objects the proper starting point 
of a nuclear many-body calculation is the translationally invari- 
ant intrinsic Hamiltonian 



Hi, 



T-T cm + V=T mt + V. 



(1) 



The use of the intrinsic kinetic energy T ml in a simple Hartree- 
Fock (HF) calculation has consequences for the validity of 
Koopmans' theorem, and thereby the interpretation of the HF 
eigenvalues as single-particle energies. A detailed analysis was 
given by Khadkikar and Kamble in Ref. |3|]and referenced re- 
peatedly over the past few decades (see e.g. 14|] or our own work 
fl5[]). However, this analysis makes explicit use of the proper- 
ties of the HF Slater determinant IT), including the assump- 
tion that it is an eigenstate of the particle-number operator with 
eigenvalue A: 

APF>=APF>. (2) 
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Naturally, this condition does not hold in a method like HFB, 
where the particle number is not conserved. We will analyze 
this case in detail in the following. 

2. Theory 

2.1. The General Case 

Since we want to deal with theories that do not conserve par- 
ticle number, we consider operators in Fock space. In this case, 
the intrinsic kinetic energy operator can be expressed either as 
a sum of one- and two-body operators, 



(3) 



or a sum of two-body operators alone, i.e., the relative kinetic 
energies of each nucleon pair: 



L y ill = _ y 



(4) 



The equality of these two expressions follows from the relation 

x2 



Ti( p '- p j) =Tj(rf +p2 j- 2prp j) 

=(^-i)2>"- 2 2>-^' (5) 



where the A resulting from the summation over the second inde- 
pendent particle i or j in the first two terms is again a Fock space 
operator measuring the total particle number. This distinction 
is inconsequential as long as one works in a Hilbert space with 
fixed particle number, because A can then be replaced by the 
corresponding eigenvalue. Naturally, one is tempted to use a 
similar replacement A — > (A) in Fock space as well, but we 
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will demonstrate in the following that this naive treatment of 
the particle-number dependence leads to discrepancies. 

Consider a many-body state \¥) without fixed particle num- 
ber. Taking the energy expectation value of the intrinsic Hamil- 
tonian with Eqs. © or © in this state, we obtain the energy 
functionals 

E (a) m = <(i - i] t) - - yA Pi ■ Pj ) + (v) (6) 



and 



£(6) m = ^Z<i(^-^) 2 > + < F >' 



(7) 



respectively. Since does not satisfy the eigenvalue equation 
I0, we have to consider the operator A directly in all expec- 
tation values, which will require a series expansion in practice. 
To this end, we note that 



A <A>1 + ^ <A)V \<A)J 



(8) 



(A) 

where we have introduced 

AA=A-(A). 

Eq. ([8]) defines a formal expansion of the energy functionals in 
powers of (Ay 1 . Applying this expansion to Zs (a) , we obtain 



(9) 



7(a) 



1 



(A) 
(TAA) 



(A)m jzj 



+ 7FrYj {Pi - p 3 hA)+0 ( {A) ~ 1 )- (10) 

(A) 2 (A) 2 m 



Denoting truncations of this series containing all terms up to 
a given order k by Ef \ the leading (LO) and next-to-leading 
order (NLO) functionals are 



Ef = (T) + (V), 



° <A> £j (A)m 



(11) 
(12) 



We note that the NLO functional is the one we would obtain by 
naively replacing A with (A) in Eq. @. 

Plugging the expansion ([8} into Eq. (|7), we have 



7 (b) 



Z< 



AA 



(A) \(A)j 



AA 



V 



2m(A) 



In this case, a naive power counting in (A) -1 breaks down be- 
cause terms at a given order of the series are enhanced by fac- 
tors of (A) as a direct consequence of applying Eq. (© to the 
expansion, i.e., 



7(b) 



(V).. 



p(b) _ p(b) 



= if + (T) - 



Z 

i<j 



({Pi-Pj) 2 ) 



2{A)m 

(T) | (TAA) 
(A) (A) 



Z 

<<j 



(Pi ■ Pj) 
(A)m 



(15) 



and 



<(/>,-/> J AA) 

2<A) 2 m 



~ (b) _ (TAA) (TAA) (T(AA) 2 ) 



z 



(A) (A) 2 
(pi ■ PjAA) 

(A) 2 m 



(A) 2 



(16) 



The term (T) that is formally of order (A) first appears in Ef\ 
a linear term in E^\ and so on. Comparing E^ with E? , we 
note that 



7(b) 



7 (a) 



(TAA) 
(A) 



(17) 



i.e., if we simply replace A with the number (A) in Eqs. © and 
(0, the functionals are no longer equivalent! 

To restore a proper power counting to the expansion of E^ b \ 
we first note that the enhanced linear term appearing in E~f 
exactly cancels the one in Eq. ( fTTl i. Likewise, an enhanced 
quadratic term in Ef cancels (T(AA) 2 )/(A) 2 in Eq. (fl6l) . and 
similar cancellations occur for all higher orders. The cancella- 
tion can be enforced explicitly if we define 



Ef = Ef + (-1) 



and applying Eq. (0 we find that 



„( T(AA) n ) 

(Ay 



p(") - p(b) 



(18) 



(19) 



i.e., the functionals E^"' and £ (/7) are identical at any given order 
of the expansion. 

2.2. Hartree-Fock-Bogoliubov Approximation 

We now apply the expansion developed in the previous sub- 
section to the HFB approximation [ 1], i.e., we assume that |*P) 
is a quasi-particle Slater determinant and introduce the density 
matrix 

pat = (T| c ; Cl m (20) 



) + (V) (13) and the pairing tensor 



(21) 



We first consider the one- plus two-body form of the intrinsic 
kinetic energy. At next-to-leading order, the functional reads 



(14) E ( ? = ll--± I -W(k\t\k')p k , k 

\ (A) I w 



1 V Pi ' P2 i 

2 Zj ^ V ~~ / Jv-" q )Pi'iP k ' k 



kk'qq' 



(A)m 



kk'qq' 



(A)m 



where t is the single-particle kinetic energy, and we assume that 
all two-body matrix elements are antisymmetrized in the fol- 
lowing. Varying Ej w.r.t. the density matrix, we obtain the 
particle-hole field 



h (a) 

"kk' 



— l — \k\ t \k') +Y(kq\ v - ^K^ 1 
(A)) jf m (A)m 

(A) 2 jzj (Afm J 



\k'q')p q >, 



where we have used 
d 1 



dpk>k (A) (A) 2 dp k , k 



(23) 



(24) 



Varying the energy with respect to the pairing tensor, we find 
that the pairing field is given by 



qq' 



Pi P2 

(A)m 



\qq')K, 



qq' ■ 



(25) 



If we start from the pure two-body form © of the intrinsic 
kinetic energy, we follow the analysis in the previous section 
and apply ( fT8l ) to obtain the NLO functional 



1 (A){rj 2p (A) 
The expectation value of the correction term is 

(TAA) = 2Y J (k\t\k')(p vk -pl k ) , 



(26) 



but due to the properties of the Bogoliubov transformation 



p — p — —kk 



(28) 



and there is an ambiguity regarding to which field the correc- 
tion term contributes when ES ^ is varied. To make contact with 
related works in the literature (see J3l and Refs. therein) as 
well as the HF treatment of Ref. Q], we split it evenly between 
the particle-hole and particle-particle channels and express the 
latter contribution in a manifestly real form: 

(TAA) = Y j (k\t\k')(p k , k -pl k ) 

kk' 

+ I J] «*| t \k')(KK*) tk + (k\ t |f> ■ (29) 

1 kk' 

We stress that the specific choice (|29] i only affects unobserv- 
able quantities like the particle-hole and pairing fields, while 
observables like the energy expectation value or the energy dif- 
ferences discussed in Sect. [3]do not depend on it at all. 
Noting that 



(Ay 



<«y> (TAA) (T) (T) 



m 



(A) 



(Pi ' Pi) 



(A) (A) 2 (A) 2 m 



(30) 



(cf. Eq.([T5ll). we find that the particle-hole field is given by 



(T) 



1 

(I) 



\t\k') -Y^{(k\t\q) Pq k' + Pk q (q\t\k')) 



(T) (T) y 



(Pi ■ Pj) 

,<A> (A) 2 j^. (A) 2 m j 
and for the pairing field, we obtain 

2 {A)m 
(A) V 

Assuming that the single-particle states satisfy 
(k\k') = So. , 



(31) 



(32) 



(33) 



we obtain the following relation for the antisymmetrized matrix 
element of the relative kinetic energy: 



-(kq\q\ 2 \k'q') 
m 

= (k\ t \k')5 qq , - (k\ t \q')8 qk , + (q\ t \q')S w 



- (q\ t \k')6 kq > (kq\ p\ ■ p 2 \k'q') ■ 

m 

(27) Plugging this into Eqs. d3~TT ) and d32b , we find that 



v(A) 



v(a) 



(34) 



(35) 
(36) 



for the NLO functionals and Ej and the specific choice 
( T29l ) for the correction term. 

The equivalence of the fields guarantees that a solution of 



the HFB equations for the functional Er will also solve the 
equations for E, and vice versa: 



[•H^ft] = [<H<» ft] =0. 



(37) 



Here, H and ft are the usual HFB Hamiltonian and generalized 
density matrices |l|], 



h-A 

-A* 



A 

-h* + A 



<R 



P 

-K* 



K 

I-P* 



and the Lagrange multiplier satisfies 



A = A {a) = A 



(38) 



(39) 



In this context, some additional remarks are in order. Since 
the expansion parameter (Ay 1 is directly linked to the vari- 
ational degrees of freedom p kk >, the derivatives d24l i generate 
(Ay 2 terms in the fields that cause global shifts in the diagonal 



matrix elements of i.e., the underlying single-particle spec- 
trum. In contrast to (A)~ 2 contributions that arise from vary- 
ing the NNLO functionals E { ^ b) , these global shifts are state- 
independent, and can be absorbed in the Lagrange multiplier A 
in a self-consistent calculation. If they are included implicitly 
in A, one cannot directly interpret A as the Fermi energy of the 
system, as it is usually done in the literature (see e.g. |7|]). 

For higher orders of the expansion, the explicit evaluation 
of the expectation values (TAA n ) and (p, ■ pjAA") occurring 
in higher-order functionals Z?„ becomes increasingly cum- 
bersome and time-consuming in practical calculations. For- 
tunately, the expansion of the energy functional converges 
rapidly, and it is sufficient to truncate the expansion of the en- 
ergy functional at the linear order in practical calculations, as 
demonstrated in Sect. [3] 



< 



10 - 




Figure 1: (Color online.) Particle number fluctuation in the tin isotopic chain. 



2.3. Hartree-Fock Limit 

Starting from the analysis of the HFB approximation in the 
previous section, we can easily take the limit of vanishing pair- 
ing correlations to obtain the Hartree-Fock limit. Since the HF 
ground state I*?) is an eigenstate of A, we immediately see that 



AA |T) = , 



(40) 



and therefore all higher-order terms in the expanded function- 
als and E (T> automatically vanish. Of course, one could 
have directly used the eigenvalue equation (f2]i and avoided the 
expansion in the first place. Moreover, this implies that the "un- 
corrected" HF functional 



- y — — + (v) 



(41) 



defined in analogy to (T5[ yields the same energies as E (a) and 
E^ b \ while we obtain the relation 

_ h (a) + 1 (f _ {tp + pf)) + <^> (42) 

for the corresponding particle-hole Hamiltonian by moving the 
"correction" terms appearing in A® to the left-hand side of Eq. 
( |35| >. Equation (l42l is exactly the relation Khadkikar and Kam- 
ble obtained by plugging ( f34b in the expression for h b From 
our analysis, we now see how this relation for the fields fol- 
lows directly from the energy functional, and that the correction 
terms in particular are derivatives of vanishing energy contribu- 
tions. 

The HFB equations (l37T i for the functionals E ia) and are 
reduced to their HF counterparts, i.e., a HF solution p obtained 
with either functional is also a valid solution for the other func- 
tional. In addition, it was demonstrated in that p also solves 
the HF equations for the uncorrected functional £ ,(i) , i.e., 



[h^,p] = [h w ,p] = [h w ,p] = 



L«0 



lb) 



(43) 



Finally, we point out one important caveat: the considera- 
tions of this subsection do not apply if the equal-filling approxi- 
mation (EFA) is employed to treat open-shell nuclei. While one 



can write down and solve HF-like equations for such a case, the 
EFA density matrix represents a statistical mixture rather than 
a genuine Slater determinant In this case, we cannot use 
the eigenvalue equation (f2]i to construct the energy functional, 
and we have to resort to the l/(A)-expansion again, treating all 
expectation values in the statistical sense. 

2.4. Particle-Number Projection 

Since nuclei are finite systems, one eventually wants to carry 
out a particle-number projection (PNP) to obtain a state with 
fixed particle number. The PNP energy functional is con- 
structed from a quasi-particle Slater determinant that is explic- 
itly projected onto the particle number A via the hermitian pro- 
jector P A 

PP A ) = P A PP> (44) 

(see e.g. |H@] for details). Since IT^) satisfies the eigenvalue 
equation ©, an expansion is not required, just as in the HF 
case. All three previously defined functionals are equivalent 
when derived from (144-b . 



p(") — p(f>) _ p\< 

^ PNP ^PNP PNP ■ 



7(b) 



(45) 



but E { p Np will in general lead to different non-observable quan- 
tities like the projected fields Iia and A^ or their individual con- 
tributions to the energy expectation value. In a variation after 
particle-number projection (VAP), the fields and densities ob- 
tained by solving the projected HFB equations |@] are associ- 
ated with an auxiliary intrinsic state without physical meaning, 
whereas the expectation values of observables are well defined, 
even if they are A-dependent. Thus, only these expectation val- 
ues of observables or derived quantities like separation energies 
(in the sense of energy differences) should be compared to ex- 
periment. 

3. Discussion & Numerical Results 

To test the proposed expansion, we have performed spheri- 
cal HFB calculations using the phenomenological Gogny D1S 
interaction ifioll . We are employing a spherical harmonic oscil- 
lator configuration space, and explicitly minimize all energies 
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Figure 2: (Color online.) Convergence of the 1 1(A) expansion: ground-state 



energy difference fiff j - 
and 120 Sn(*). 



pto/f") 



as a function of the order r'for ls O (•), 64 Ni(H), 



with respect to the oscillator length aho; more details can be 
found in Ref. JfJ. Unless noted otherwise, we include 13 major 
oscillator shells in our calculations, which leads to a satisfactory 
convergence in the considered cases. Odd nuclei are treated in 
a self-consistent blocking method yj], where the odd nucleon is 
distributed evenly over all magnetic substates of a given y'-shell 
according to the equal-filling approximation (see e.g. HI). 

Let us first consider the convergence of the l/(A)-expansion. 
Formally, an A-quasi-particle HFB state can contain states with 
sharp particle numbers from 0, . . . , 2A. In the extreme cases, 
this would mean that the operator AA/(A) in Eq. ([8]) could 
acquire the operator norm 1 on the space of HFB wavefunc- 
tions, causing the breakdown of the series expansion. In the 
nuclear many-body problem, this breakdown could only oc- 
cur before the first major shell is fully occupied, and for these 
very light nuclei the use of mean-field methods is questionable 
in the first place. At the major shell closure itself, the HFB 
wavefunction collapses onto the HF solution, and there is no 
need for an expansion. As we progress along the nuclear chart, 
we find that particle-number fluctuations only occur within a 
given major shell, as shown exemplary for the tin isotopes with 
N = 50, . . . , 82 in Fig. Q] This implies that the operator norm 
of AA/(A) remains below 1 in practical applications, and guar- 
antees the convergence of the series. 

To test the rate of the convergence, we display the quantity 

E (alb) _ E (a/b) for open . she ii nucle j 18q 64^ ^ 120 Sn j n 

Fig. |2] We have picked these specific nuclei for their large val- 
ues of (AA 2 ) in the respective isotopic chains. We find that is 
essentially sufficient to include the linear terms in 1/(A) in the 
energy functional. Beyond the linear order, the largest varia- 
tions occur for 18 0, where the successive inclusion of terms up 
to third order causes changes of 100 - 200 keV (see the inset 
of Fig. 0. As expected, the effect of higher orders rapidly di- 
minishes with increasing masses, and amounts to a few keV for 
the nickel and tin isotopes. For this reason, we will consider all 
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Figure 3: (Color online.) Theoretical three-point binding energy differences in 
the tin isotopic chain from £w" (#) and £W (■) t compared to experimental 
values (+)(jj]. 



functionals at linear (i.e. next-to-leading) order, and drop the 
subscripts in the following. 

Information about the effect of the intrinsic kinetic energy on 
the nuclear pairing correlations can be extracted by considering 
the three-point binding energy difference formula 

A (3) (A0 = (.-lf-r (E(N + 1) - 2E(N) + E(N - 1)) (46) 

along an isotopic chain (or equivalently, an isotonic chain with 
N replaced by Z). By calculating the ground states of odd nuclei 
self-consistently, we avoid the complications arising in pertur- 
bative analyses, as discussed in detail in yjl I13H. This is par- 
ticularly relevant for A-dependent Hamiltonians |6i], which add 
another layer of complexity to the perturbative approach. 

In Fig. [3j we compare the theoretical A (3) for the properly 
constructed functional to the uncorrected functional E (b> ( IT5l . 
The latter are typically lower than the values for the proper 
functional by 200 - 400 keV, with the exceptions occurring 
around the W5/2 and W3/2 sub-shell closures at N = 56 and 
N = 70, respectively. The gaps between the relevant sub-shell 
and the next available one are notably larger for E®> than for 
E ( - a '° ) , leading to more pronounced effects when the odd nu- 
cleon is added 

In Sect. 12.31 we have pointed out that a spherical HF treat- 
ment of open-shell nuclei will suffer from the same problems as 
the HFB method due to the use of the equal-filling approxima- 
tion. This is explicitly demonstrated for the tin isotopic chain 
in Fig. |4] where we have used the density matrix p obtained by 
minimizing the functional E^ F to calculate the energies E^ HF 
and E^ p . At sub-shell closures, we have a genuine HF prob- 
lem, and in this case all three functionals are equivalent and p 
is a solution for each set of HF equations [3]. For open-shell 
nuclei, in contrast, the binding energies obtained with E^ F are 
reduced by several hundred ke V, and the difference is given by 



<^AA) 
(A) 



= — tr (p-p 2 ) 
(A) ^ H) 



(47) 



to numerical accuracy. 
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Figure 4: (Color online.) Ground-state energy differences for the tin isotopes in 
HF+EFA: £™ - E<g> p (•) and E®> p - E$ p (■). Calculations were done with 
1 1 major shells. 
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Figure 5: (Color online.) Top: Ground-state energy differences E^^ p - EpL p 
of the tin isotopes from VAP calculations (note the scale in eV). Bottom: Total 
pairing energies of the tin isotopes from VAP calculations with Ep% p (•) and 
Epnp (^)- All calculations were done with 1 1 major shells. 



of the pairing energies in the bottom half of Fig. illustrates 
that the individual particle-hole and particle-particle contribu- 
tions to the total energy expectation value may well be different 
in both calculations. 

4. Conclusions 

We have presented a detailed analysis of the A-dependent 
intrinsic kinetic energy operator in theories with and without 
particle-number conservation, in particular the mean-field type 
Hartree-Fock-Bogoliubov method and its number-conserving 
extension via particle-number projection. We have shown that a 
naive treatment where the number operator A is replaced by its 
expectation value (A) causes discrepancies in expectation val- 
ues obtained with the otherwise equivalent operator forms of 
T mt . We have developed a systematic expansion to fix this prob- 
lem, but we emphasize that this expansion does not restore, nor 
is it intended to restore, either the particle number or transla- 
tional symmetries of the nucleus. 

Our discussion provides an a posteriori justification for us- 
ing the one- plus two-body form of the intrinsic kinetic energy 
since it is automatically consistent with the power counting of 
the developed expansion. As a byproduct, we also clarify how 
differences and ambiguities in non-observable quantities which 
had been discussed in the context of the HF approximation J3t] 
arise systematically in the presented framework. While we have 
discussed the specific case of the intrinsic kinetic energy oper- 
ator in the present article, we point out that the same treatment 
should be applied to all A-dependent observables. 
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Finally, we compare the tin ground-state energies obtained References 



with Ep Np to those of E^p' in calculations where the variation 
is carried out after particle-number projection (VAP). We are 
aware that the VAP method is highly problematic for density- 
dependent interactions like Gogny D1S (see e.g. IU4I1 and ref- 
erences therein). Since we do not aim for a comparison with 
experiment in the present discussion, we assume that the po- 
tential energies of both calculations will be equally affected by 
any spurious behavior due to the density-dependence, hence 
any differences in the VAP energy expectation values are due 
to the different forms of the kinetic energy. Looking at the top 
half of Fig. [5] and noting the eV scale, we find no such dif- 
ferences: the tin ground-state energies obtained in VAP calcu- 
lations with either the properly constructed functional EpLp or 
the naive, uncorrected functional Ep Np are identical, confirming 
the discussion of Sect. 12.41 At the same time, the comparison 
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